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I  INTRODUCTION 


The  disturbed,  nighttime  equatorial  ionosphere  has  been  found  to 
contain  two  unique  features,  "plasma  bubbles"  and  radar  backscatter 
"plumes."  Plasma  bubbles  are  localized  depletions  in  F- layer  electron 
density  which  can  be  as  much  as  three  orders  of  magnitude  below  the 
ambient  level  (Hanson  and  Sanatani,  1973;  McClure  et  al.,  1977).  Plasma 
bubbles  are  believed  to  extend  throughout  entire  magnetic  flux  tubes 
v/ith  dimensions  transverse  to  the  geomagnetic  field  of  up  to  a  few 
hundred  kilometers.  Radar  backscatter  plumes  are  produced  by  meter- 
scale  field-aligned  irregularities  (FAI)  and  extend  from  the  bottomside 
into  the  topside  of  the  F  layer.  Because  plasma  bubbles  and  radar  plumes 
appear  to  possess  similar  characteristics,  researchers  have  speculated 
that  they  are  directly  associated,  if  not  spatially  coincident  (e.g., 
Woodman  and  La  Hoz,  1976). 

Although  there  have  been  several  coordinated  experiments  (involving 
in-situ  measurements  of  plasma  depletions  by  rockets)  whose  results  have 
suggested  that  a  close  relationship  does  exist  (Kelley  et  al.,  1976; 

Morse  et  al.,  1977),  the  results  have  not  been  conclusive.  Perhaps  the 
most  convincing  evidence  published  thus  far  are  those  of  Tsunoda  and 
Towle  (1979).  They  showed  that  large-scale  (200  km  in  width)  depletions 
in  total  electron  content  (TEC)  were  longitudinally  coincident  with 
radar  backscatter  plumes.  However,  because  TEC  in  a  range- integrated 
measurement,  it  was  not  possible  to  determine  whether  the  regions  of 
plasma  depletion  that  produced  the  TEC  depletions  were  spatially  coinci¬ 
dent  with  the  backscatter  plumes,  or  whether  they  were  displaced  in 
altitude  from  one  another. 

In  this  paper,  we  present  spatial  maps  of  plasma  bubbles  obtained 
with  an  incoherent-scatter  (IS)  radar.  (Similar  results  have  also  been 
reported  by  Towle  [1979]).  Using  ALTAIR,  an  IS  radar  located  in  the 
Kwajalein  Atoll,  Marshall  Islands  (Tsunoda  et  al.,  1979),  both  plasma 
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bubbles  and  radar  plumes  were  spatially  mapped,  sequentially  in  time. 

A  comparison  of  the  locations  of  plasma  depletion  and  radar  backscatter 
indicate  that  they  are  directly  related. 


II  THE  EXPERIMENT 


The  use  of  ALTAIR,  a  fully- steerable  incoherent- scatter  (IS)  radar, 
for  the  spatial  mapping  of  equatorial  field-aligned  irregularities  (FAI) 
and  the  electron-density  distribution  of  the  background  F  layer  was  de¬ 
scribed  by  Tsunoda  et  al.  (1979).  Briefly,  ALTAIR  is  located  in  the 
Kwajalein  Atoll,  Marshall  Islands  at  a  magnetic  dip  latitude  of  4.3°N. 
The  radar  was  operated  for  this  experiment  at  155.5  MHz,  transmitting  a 
30-ps  pulse  through  a  46-m  paraboloid  antenna  having  a  2.8°  beamwidth. 

In  order  to  map  plasma  depletions  using  incoherent  scatter,  we  must 
be  able  to  suppress  the  contaminating  effects  of  spread-F  backscatter. 
Suppression  of  intense  spread-F  backscatter  is  accomplished  by  taking 
advantage  of  the  field- aligned  nature  of  the  irregularities  and  the 
steerability  of  the  ALTAIR  antenna.  Because  of  the  strong  magnetic 
aspect  dependence  of  spread-F  backscatter,  IS  measurements  can  be  made 
with  excellent  suppression  of  spread-F  backscatter  by  positioning  the 
radar  beam  several  degrees  away  from  perpendicular  aspect  with  the  geo¬ 
magnetic  field. 

A  second  consideration  is  the  capability  of  the  radar  to  resolve  a 
plasma  bubble.  The  smallest  plasma  bubble  that  can  be  resolved  with 
ALTAIR  is  dependent  on  the  range  to  the  depletion.  At  400  km,  a  typical 
distance  from  radar  to  a  plasma  bubble,  the  dimensions  of  the  radar 
scattering  volume  are  4.5  X  20  x  20  km.  Therefore,  ALTAIR  is  capable  of 
spatially  mapping  a  plasma  bubble  that  is  larger  than  4.5  km  in  range 
and  20  km  in  transverse  distance.  It  is  worth  noting  that  the  altitude 
extent  of  plasma  depletions  observed  to  date  by  rockets  is  on  the  order 
of  10  km,  or  less  (Kelley  et  al.,  1976;  Morse  et  al.,  1977).  On  the 
other  hand,  the  horizontal  dimensions  of  plasma  bubbles  appear  to  extend 
up  to  200  km  (McClure  et  al.,  1977), 

The  experiment,  conducted  on  18  August  1978,  consisted  of  operating 
ALTAIR  in  two  discrete-scan  modes,  with  one  following  immediately  after 


the  other.  The  scan  geometry  for  the  two  modes  is  shown  in  Figure  1, 

The  coordinates  are  polar  with  radar  azimuth  measured  from  true  north 
(top  of  the  figure)  in  a  clockwise  direction.  The  first  mode  consisted 
of  making  a  discrete,  25- position  scan  from  west  to  east  over  a  72° 
angular  sector,  maintaining  the  radar  beam  perpendicular  to  the  geomagnetic 
field  at  F-region  altitudes.  The  odd-numbered  beam  positions  for  this 
scan  mode  are  shown  in  Figure  1.  The  scan  is  seen  to  intersect  the 
magnetic  meridian  at  approximately  9°  north  of  vertical.  The  radar  beam 
was  held  fixed  for  about  15  seconds  in  each  position,  thus  completing 
the  scan  in  six  minutes.  With  this  type  of  scan,  we  were  able  to  map  the 
spatial  distribution  of  radar  backscatter  plumes. 

Upon  completion  of  the  first  scan  mode,  the  location  of  radar  plumes 
was  determined  from  real-time  range- time- intensity  (RTI)  displays.  Five 
off-perpendicular  beam  positions  were  then  selected  for  the  second  scan 
mode  which  corresponded  (when  extrapolated  along  the  magnetic  field  lines) 
to  the  vicinity  of  maximum  backscatter  strength  within  a  radar  plume. 

The  five  beam  positions  were  selected  from  25  preselected,  off-perpen¬ 
dicular  beam  positions  that  corresponded  to  those  in  the  first  scan  mode, 
as  shown  in  Figure  1.  Incoherent- scatter  measurements  were  then  recorded 
for  one  minute  at  each  of  those  beam  positions. 

The  off-perpendicular  beam  positions  in  Figure  1  (primed  numbers) 
correspond  to  the  beam  positions  in  the  first  scan  mode  labeled  with  the 
same,  but  unprimed  numbers.  Note  that  the  off-perpendicular  beam  positions 
are  offset  southward  from  the  unprimed  beam  positions  by  about  13°.  With 
this  displacement,  the  first  sidelobe  (estimated  to  be  approximately  6° 
away  from  boresight  of  the  main  beam)  is  kept  well  away  from  the  contour 
of  exact  orthogonality  with  the  magnetic  field  lines.  Consequently,  we 
believe  that  suppression  of  backscatter  from  spread-F  irregularities  is 
well  below  -36  dB,  the  (two-way)  gain  of  the  first  sidelobe  relative  to 
the  main  lobe.  In  fact,  the  suppression  is  probably  closer  to  -50  dB, 
as  indicated  by  the  data  presented  in  the  following  section. 

In  order  to  assure  this  level  of  suppression,  we  sacrificed  near- 
spatial  ly-coincident  measurements  and  introduced  the  assumption  that 
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FIGURE  1  RADAR  SCAN  GEOMETRY  SHOWING  THE  BEAM  POSITIONS  OF  THE  SCAN  MADE  PERPENDICULAR 
TO  THE  MAGNETIC  FIELD  AND  THE  CORRESPONDING  OFF-PERPENDICULAR  BEAM  POSITIONS 


* 


both  1-m  FAI  features  and  plasma  depletions  map  along  magnetic  field 
lines.  This  assumption  is  considered  to  be  very  reasonable  since  13° 
angular  separation  at  400-km  range  corresponds  to  a  spatial  separation 
of  90  km,  or  less  than  a  degree  of  latitude.  Equatorial  spread-F  pheno¬ 
mena  are  known  to  occur  over  ±20°  in  latitude  and  are  believed  to  involve 
entire  magnetic  flux  tubes.  Morse  et  al.  (1977)  found  that  ionospheric 
features  tended  to  be  similar  over  several  degrees  of  latitude.  Un¬ 
published  ALTAIR  IS  data,  collected  during  elevation  scans  in  the  magnetic 
meridian,  also  support  this  hypothesis. 
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Ill  RESULTS 


The  spatial  distribution  of  backscatter  from  1-m  FAI  obtained  from 
the  first  scan  mode  is  shown  in  Figure  2.  The  contours  of  constant 
backscatter  strength  (corrected  for  range  dependence)  which  are  plotted 
in  10-dB  steps,  are  referenced  to  the  IS  signal  strength  associated  with 
an  electron  density  of  10^  el/criP.  In  other  words,  the  0-dB  contour 
corresponds  in  signal  strength  to  incoherent  backscatter  that  would  be 

/•  O 

obtained  from  an  electron  density  of  10°  el/cm  .  The  10-dB  contour 
corresponds  to  an  equivalent  electron  density  of  10^  el/ctip,  and  so  on. 

Two  backscatter  plumes  are  seen  in  Figure  2,  separated  by  a  region 
of  no  backscatter.  The  larger  plume  in  the  west  half  of  the  contour  map 
is  seen  to  be  connected  to  backscatter  from  the  bottomside  F  layer  and 
extends  upward  (and  westward)  to  altitudes  near  700  km.  Both  radar 
plumes  are  typical  of  those  seen  with  ALTAIR  on  very  disturbed  nights 
(e.g.,  Tsunoda  and  Towle,  1979)  and  similar  to  those  observed  at  50  MHz 
with  the  Jicamarca  IS  radar  (Woodman  and  La  Hoz,  1976). 

The  region  of  strongest  backscatter  (greater  than  30  dB)  occurred 
in  beam  positions  5  through  8.  Having  examined  the  RTI  display  (corres¬ 
ponding  to  Figure  2)  in  real  time,  we  estimated  the  locations  of  strongest 
backscatter  and  selected  of f- perpendicular  beam  positions  6 '  through  107 
for  IS  measurements. 

The  electron  density  profiles  obtained  from  the  five  of f- perpendicular 
beam  positions  are  shown  in  Figure  3.  Each  profile  is  labeled  with  its 
time  of  measurement,  and  shaded  in  regions  where  the  electron  density  is 
greater  than  10“*  el/cnP.  Extreme  care  should  be  taken  in  the  interpretation 
of  the  structure  seen  in  the  profiles,  particularly  since  the  IS  data 
were  noisy  above  600  km.  Because  of  poor  signal- to- noise  ratio,  the 
variations  in  electron  density  above  600  km  should  not  be  considered  to 
be  real.  Together  with  the  profiles  are  shown  vertical  black  bars  which 
represent  the  altitude  extent  over  which  spread-F  backscatter  was  observed. 


11 


ANTENNA  POSITIONS 


12 


13 


Examining  the  electron  density  profiles,  we  find  that  all  of  the 
profiles  are  characterized  by  a  major  depletion  located  near  450  km 
altitude,  very  close  to  the  peak  of  the  F  layer.  The  plasma  depletion 
extends  over  100  km  in  altitude  with  a  percentage  depletion  that,  in 
some  cases,  approaches  two  orders  of  magnitude.  All  of  the  depletions 
also  have  similar  internal  structure,  i.e.,  two  minima  in  electron  density, 
suggesting  that  perhaps  there  are  actually  two  depletions  contiguous  in 
altitude . 

The  black  vertical  bars  shown  together  with  each  electron  density 
profile  actually  represent  our  interpolated  estimates  of  where  spread-F 
backscatter  must  have  occurred  at  the  times  the  electron  density  measure¬ 
ments  were  made.  Interpolation  is  necessary  because  of  an  eastward, 
bulk-plasma  drift  that  prevails  in  the  nighttime  equatorial  ionosphere 
(Woodman,  1972).  Without  actual  plasma  drift  measurements,  we  used  a 
nominal  100  m/s  eastward  drift  speed  as  our  best  estimate.  Then  using 
the  horizontal  spacing  of  the  electron  density  profiles  in  Figure  3  as 
the  relative  spacing  between  beam  positions,  we  interpolated  the  locations 
of  spread-F  backscatter  measured  during  the  first  scan,  made  perpendicular 
to  the  magnetic  field. 

The  vertical  lines  labeled  by  numbers  in  Figure  3  show  the  relative 
locations  of  backscatter  features  measured  at  those  beam  positions,  at 
times  corresponding  to  the  off-perpendicular  measurements.  For  example, 
the  backscatter  altitudes  corresponding  to  the  electron  density  profile 
taken  at  1155:30  UT  were  estimated  from  the  backscatter  in  beam  positions 
4  (not  shown)  and  5.  The  backscatter  altitudes  associated  with  the  pro¬ 
file  taken  at  1156:45  UT  were  interpolated  from  beam  positions  5  and  6 
that  bracket  the  profile  in  Figure  3.  The  backscatter  altitudes  asso¬ 
ciated  with  the  remaining  three  electron  density  profiles  were  estimated 
from  beam  positions  6,  7,  and  8,  respectively.  Except  for  the  backscatter 
corresponding  to  the  last  electron  density  profile,  the  interpolation 
was  not  critical  in  placing  the  backscatter  regions. 

In  all  cases,  we  see  that  backscatter  is  associated  with  (1)  the 
bottomside  F  layer,  and  (2)  the  plasma-depleted  region.  The  only  other 
backscatter  that  was  observed  were  those  found  in  the  topside  layer  at 
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1200  UT.  As  shown,  the  spatial  coincidence  of  radar  backscatter  and 
plasma  depletions  is  excellent.  The  altitude  extent  of  backscatter  is 
seen  to  be  comparable  to  the  altitude  extent  of  the  plasma  depletion. 

The  fact  that  they  both  occur  at  similar  altitudes  and  have  comparable 
thicknesses  (in  altitude)  is  perhaps  the  most  convincing  evidence  that 
radar  backscatter  and  plasma  bubbles  are  directly  related. 

If  we  compare  the  altitudes  of  the  plasma  depletions  found  in  each 
profile,  we  find  that  the  depletions  occur  at  lower  altitudes  as  we  look 
further  eastward.  Since  the  five  beam  positions  are  contiguous,  the 
systematic  variation  in  depletion  altitude  can  be  interpreted  as  evidence 
of  a  single,  slab- like  depletion  in  electron  density  that  is  tilted. 
Ignoring  any  vertical  transport  that  might  be  associated  with  the  plasma 
depletion,  such  as  an  upward  bubble  rise- velocity  (e.g.,  Ossakow  and 
Chaturvedi,  1978),  or  a  bulk  downward  movement  of  the  entire  F  layer,  we 
find  that  the  slab- like  depletion  is  tilted  at  an  apparent  angle  of  60° 
west  of  vertical. 

Referring  back  to  Figure  2,  we  note  that  the  tilt  found  in  the 
plasma  depletion  is  in  the  same  sense  as  the  tilt  found  in  the  20  dB 
contour  of  corresponding  backscatter  region.  This  comparison  is  reason¬ 
able  because  the  location  and  altitude  extent  of  the  20  dB  contour  are 
comparable  to  that  of  the  plasma  depletion.  However,  a  quantitative 
comparison  of  tilt  angles  is  difficult.  (For  example,  the  tilt  of  the 
30  dB  contour  is  completely  different  from  that  of  the  20  dB  contour.) 

The  exact  spatial  relationship  between  the  two  phenomena  is  less 
clear.  If  we  extrapolate  from  the  region  of  spread-F  backscatter  along 
a  magnetic  flux  tube  to  the  location  of  IS  measurement,  we  must  shift 
the  backscatter  region  upward  in  altitude  by  about  15  km.  We  must  also 
account  for  any  vertical  transport  that  might  have  occurred  between  the 
time  of  spread-F  backscatter  measurements  and  the  time  of  IS  measurements. 

For  a  more  detailed  comparison  between  plasma  bubbles  and  spread-F 
backscatter,  we  have  plotted  the  backscatter  profiles  from  beam 
positions  5  through  8  together  with  their  corresponding  electron  density 
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profiles,  in  Figure  4.  The  backscatter  from  plume  FAI  is  seen  to  have 
two  maxima  in  all  four  profiles,  suggesting  a  relationship  to  the 
internal  structure  of  the  plasma  depletions. 

As  plotted,  the  backscatter  maxima  tend  to  occur  at  altitudes 
corresponding  either  to  the  electron  density  minima  or  the  upper  walls 
of  the  plasma  depletions.  If  we  allow  for  the  15-km  altitude  correction 
based  on  extrapolation  along  magnetic  field  lines,  the  backscatter  maxima 
would  appear  to  be  associated  with  the  upper  walls  of  the  plasma  de¬ 
pletions.  However,  without  knowledge  of  the  vertical  transport  velocity, 
it  is  difficult  to  conclude  what  the  exact  spatial  relationship  is  be¬ 
tween  radar  backscatter  maxima  and  plasma  depletions. 

Although  vertical  transport,  in  particular,  that  due  to  an  upward 
bubble  velocity,  could  drastically  change  the  derived  spatial  relation¬ 
ship  between  radar  backscatter  and  the  plasma  depletion,  the  backscatter 
profiles  in  Figure  4  can  be  used  to  argue  that  the  backscatter  cannot  be 
grossly  displaced  with  respect  to  the  plasma  depletions.  Tsunoda  et  al . 
(1979)  showed  that  IS  measurements  can  be  made  even  when  the  radar  beam 
is  directed  perpendicular  to  the  magnetic  field,  provided  that  spread-F 
backscatter  is  absent.  Therefore,  if  a  plasma  depletion  existed,  without 
spatially-coincident  spread-F  backscatter,  the  depletion  should  be  evident 
in  the  backscatter  profiles  in  Figure  4.  It  is  clear  from  examining  these 
profiles  that  the  depletions  are  not  significantly  displaced  with  respect 
to  the  backscatter  region.  From  this  observation,  we  conclude  that  the 
upward  bubble  velocity  must  have  been  negligibly  small  and  did  not  distort 
the  spatial  relationship  between  spread-F  backscatter  and  plasma  depletions. 

To  lend  further  credence  to  the  results  just  described,  let  us 
examine  other  electron  density  profiles  obtained  with  the  radar  beam 
fixed  at  10°  south  (geomagnetic)  of  vertical.  In  Figure  5,  we  show  two 
electron  density  profiles  taken  shortly  after  the  previous  data  set.  On 
the  basis  of  the  assumed  100  m/s  eastward  drift  speed,  these  profiles 
would  correspond  to  features  seen  in  Figure  2  that  are  located  about 
100  km  west  of  ALTAIR.  Both  profiles  show  a  plasma  depletion  at  altitudes 
below  400  km,  altitudes  that  seem  to  correspond  to  the  leading  edge  of 
the  20-dB  contour  in  Figure  2.  Above  the  depletions,  the  profiles  extend 
to  altitudes  above  600  km. 


ALTITUDE  —  km  ALTITUDE 


Lxwm  IS  PROFILE 


ELECTRON  DENSITY  —  el/cm3  ELECTRON  DENSITY  —  el/cm3 


FIGURE  4  A  COMPARISON  OF  PLUME  BACKSCATTER  PROFILES  (solid  curves) 

AND  CORRESPONDING  ELECTRON  DENSITY  PROFILES  (dashed  curves) 
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ALTITUDE 


FIGURE 


ELECTRON  DENSITY  —  «l/em3 


5  ELECTRON-DENSITY  PROFILES  OBTAINED  ON  18  AUGUST  1978 

WITH  RADAR  BEAM  FIXED  AT  10°  SOUTH  (geomagnetic)  OF  VERTICAL 
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The  two  profiles  in  Figure  5  can  be  compared  to  another  pair  of 
electron-density  profiles,  shown  in  Figure  6,  taken  at  the  same  time  on 
the  following  night  when  spread-F  activity  was  much  weaker.  The  profiles 
in  Figure  6  display  a  peak  electron  density  that  appears  to  be  somewhat 
higher  than  those  in  Figure  5,  although  the  depletions  (which  occurred 
very  close  to  the  altitude  of  the  F- layer  peak)  obscure  the  true  peak 
electron  density.  The  layer  thickness  of  the  F  layer  on  both  nights 
was  comparable,  with  the  layer  extending  from  250  to  650  km  in  altitude. 
The  only  noticeable  difference  in  the  background  profiles  is  the  slope 
of  the  topside  of  the  F  layer,  being  steeper  on  the  quieter  night. 
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E 


FIGURE 


6  ELECTRON-DENSITY  PROFILES  OBTAINED  ON  19  AUGUST  1978 

WITH  RADAR  BEAM  FIXED  AT  10°  SOUTH  (geomagnetic)  OF  VERTICAL 


IV  DISCUSSION  AND  CONCLUSIONS 

The  results  presented  in  this  paper  are  some  of  the  first  on  the 
two-dimensional  spatial  mapping  of  altitude-extended  plasma  bubbles. 
(Similar  results  have  been  reported  by  Towle  [1979].)  In-situ  satellite 
measurements  (Hanson  and  Sanatani,  1973;  McClure  et  al.,  1977)  have  all 
been  longitudinal  cuts  through  plasma  bubbles.  Plasma  bubbles  penetrated 
by  rockets  (Kelley  et  al.,  1976;  Morse  et  al.,  1977)  have  all  been  no 
more  than  10-20  km  in  altitude  extent  with  depletions  of  approximately  50%. 

The  plasma  bubbles  presented  in  this  paper  extended  over  100  km  in 
both  altitude  and  longitude,  and  had  percentage  plasma  depletions  ap¬ 
proaching  99%.  The  factor  that  probably  limits  the  highest  percentage 
depletion  that  can  be  measured  is  the  antenna  sidelobes.  Our  estimate 
of  -50-dB  effective  sidelobe  gain  appears  to  be  confirmed  by  the  data  in 
Figure  4.  The  maximum  separation  between  minima  in  plasma  depletions 
and  maxima  in  backscatter  strength  is  no  more  than  50  dB. 

The  comparison  of  1-m  FAI  distribution  with  the  electron  density 
profiles  associated  with  plasma  bubbles  revealed  a  direct  spatial  rela¬ 
tionship.  Backscatter  was  found  (1)  to  occur  at  the  same  altitude  and 
(2)  to  have  similar  altitude  extent  as  plasma  depleted  regions.  Two 
maxima  in  backscatter  strength  were  found  to  correspond  with  two  minima 
in  electron  density.  The  only  uncertainty  is  whether  the  backscatter 
is  associated  with  the  walls  of  the  plasma  bubbles  or  whether  it  occurs 
throughout  the  entire  bubble  volume.  Because  both  backscatter  and  plasma 
depletions  have  comparable  altitude  extents,  it  is  unlikely  that  radar 
backscatter  is  always  confined  to  steep  gradients  associated  with  the 
walls  of  the  plasma  bubbles.  Plasma  structure  must  also  exist  within 
the  bubble  volume  that  is  capable  of  producing  significant  backscatter. 

The  electron  density  profiles  shown  in  Figure  5  suggest  that  a 
plasma  depletion  is  not  only  associated  with  the  "head"  of  the  radar 
backscatter  plume  (i.e.,  the  most  intense  portion  of  the  backscatter 
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region  that  usually  appears  at  high  altitudes)  but  also  with  the  "neck" 
of  the  radar  plume  (i.e.,  the  portion  that  connects  the  plume  head  to 
the  bottomside  of  the  F  layer).  This  conclusion  supports  the  correct¬ 
ness  of  the  collisional  Rayleigh-Tay lor  instability  model  as  the  principal 
source  mechanism  of  plasma  bubbles  (e.g.,  Ossakow  and  Chaturvedi,  1978). 
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